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Abstract 

In 2006, a simple extension of the Standard Model was proposed in which neutrinos 
obtain radiative Major ana masses at one- loop level from their couplings with dark 
matter, hence the term "scotogenic," from the Greek "scotos" meaning darkness. Here 
an analogous mechanism for Dirac neutrino masses is discussed in a minimal model. In 
different ranges of the parameter space, various candidates for dark matter are possible. 
In particular, the lightest Dirac fermion which appears in the loop diagram generating 
neutrino mass can be a viable dark matter candidate. Such a possibility does not exist 
for the Majorana case. Realistic neutrino mixing in the context of A4 is discussed. A 
possible supersymmetric extension is also briefly discussed. 



Dirac neutrino masses have not received much attention in the literature mainly because 
of two reasons: (1) In the Standard Model (SM) of particle interactions, there are left-handed 
lepton doublets (u, I) l and right-handed charged-lepton singlets Ir but no v R because it 
transforms trivially under the SU(3)c x SU{2) L x U(l)y gauge symmetry and there is no 
need for its existence. If it is added in by hand, the neutrino can then obtain a Dirac mass mr, 
in the same way as all the other fermions (quarks and charged leptons), i.e. from the vacuum 
expectation value of the scalar Higgs doublet of the Standard Model. However, since vr is a 
neutral singlet, there is no symmetry which prevents it from having a large Majorana mass 
M. As a result, obtains an effective small Majorana mass from the seesaw mechanism [Tj, 
i.e. m u ~ —m 2 D /M. (2) If a symmetry is imposed in such a way that the lepton number is 
conserved, the Majorana mass term for vr will be forbidden. In that case, because neutrino 
masses are known to be of order 1 eV or less, the corresponding Yukawa couplings must 
be of order 10~ n or smaller. Such a small value is considered by many to be intrinsically 
unacceptable. 

Nevertheless, up to now, there is not any indisputable evidence for the Majorana nature of 
the neutrinos from the searches for the neutrinoless double beta decay. Thus, the possibility 
of Dirac neutrino masses cannot be discounted. To overcome the above theoretical objections, 
it is proposed in this paper that neutrinos are Dirac fermions, with two important properties. 
(1) They are protected from becoming Majorana fermions by a U(1)b-l global or gauge 
symmetry. (2) They are protected from having a tree-level mass by a Z 2 symmetry which 
is identifiable with that of dark matter, as well as another Z 2 symmetry which sets them 
apart from other Dirac fermions. The latter symmetry is broken explicitly by soft terms. 
It may also be replaced by supersymmetry, but that would require a much larger Higgs 
content. As a result, neutrinos acquire one- loop radiative masses through their couplings 
with dark matter, hence the term "scotogenic," from the Greek "scotos" meaning darkness. 
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These Dirac neutrino masses can be highly suppressed in the same way that the usual seesaw 
Majorana neutrino masses are highly suppressed. Their smallness can be also explained by 
the smallness of the soft terms breaking the Z 2 symmetry. 

In 2006, it was proposed [2j[3] that neutrinos are massive only because of their couplings 
with dark matter. This idea connects two of the most important issues in the particle physics 
and astrophysics. The idea was easily implemented [2] in a simple extension of the Standard 
Model by adding a second scalar doublet (rj + , 7]°) and three neutral singlet fermions Ni which 
are odd under an extra exactly conserved discrete Z 2 symmetry |3], in analogy to the R parity 
in supersymmetry. As a result, either t)r = y/2Re(r]°) or the lightest N may be considered 
a candidate for dark matter. In particular 7] has been called the "inert" scalar doublet in a 
model proposed [5] after Ref. [2] and studied by many authors since then [6]. Variations of the 
original idea also abound and have become an active area of research [7J |9l [101 HH [121 03] • 

In almost all previous such applications, neutrino masses have always been assumed to 
be of Majorana type. Suppose they are exactly Dirac. Is the connection between neutrino 
mass and dark matter still possible? If so, what are the necessary theoretical ingredients 
for it to happen, and what are the phenomenological consequences? In [14] . using scalar 
singlets, a radiative Dirac neutrino mass is obtained; however, in this mechanism, the dark 
matter fields do not propagate in the loop. Employing an idea similar to that proposed 
in Ref. [2], Ref [15] suggests a model both for a dark matter candidate and generation of 
radiative Dirac neutrino mass. As indicated below, this model shares some features with the 
model introduced in the present paper. In [13], a model is introduced in which neutrinos 
obtain a Dirac mass via a one-loop diagram similar to that in [15] and a Majorana mass 
via two-loop diagrams after spontaneous breaking of the lepton number symmetry. In our 
model described below, the neutrino mass is purely of the Dirac type. 

Consider first the imposition of a conserved additive lepton number to protect the neu- 
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trino mass from becoming Majorana. We choose to do so by extending the Standard Model 
to include B — L as either a global or gauged U(l) symmetry. The latter has long been 
known to be a well-motivated anomaly-free extension which requires the existence of three 
singlet right-handed neutrinos. Of course, in breaking the gauged U(1)b-l, we have to be 
sure that the global U(1)b-l symmetry of the sector relevant to the present study remains 
intact. This can be done easily by a scalar field transforming under U(1)b-l but not cou- 
pling to other fields with nonzero B — L. The second step is to forbid a tree-level Dirac 
neutrino mass m„, and yet allow a tree- level charged- lepton mass mi. To do this, the sim- 
plest way is to impose a Z% symmetry such that v c is odd but all other fermions are even. 
There is therefore no connection between v and v c at the tree level. To make them connect 
in one loop, new particles are postulated which are odd under an exactly conserved Z\ , 
and the previous Z% is allowed to be broken by soft terms. Another way is to make the 
model supersymmetric as well so that mi comes from $x = (fa fa) but m v is forbidden to 
couple to $2 = (<^2~)02) wn i cn i s assumed odd under Z^\ In either case, we need to add 
heavy neutral singlet Dirac fermions (Ni,Nf) of odd Z^ transforming under U(1)b-l and 
a neutral singlet scalar x° of odd Z^ which is trivial under U(1)b-l- 

First, let us consider the minimal non-supersymmetric model. It is a simple extension 
of the Standard Model in the same spirit of Ref. [2j. Its particle content is listed in Table 
1. In addition to the usual particles of the Standard Model, we have added three copies of 
the Weyl spinors z/ c , three copies of the Dirac spinor pairs (N, N c ), one extra scalar doublet 
V = (v + >V°) an d one real scalar x ■ The B — L symmetry prevents N, N c as well as v c 
from having a Majorana mass. Note that Z^ is broken softly by the trilinear term Ax&m 
whereas Z^ remains unbroken. ($ = (fa , <f)°) is the standard model Higgs doublet.) The 
one- loop Dirac neutrino mass is thus generated, as shown in Fig. 1. 

Note that x° is essential here for a scotogenic Dirac neutrino mass, whereas the scalar 
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Table 1: Assignments of the particles of the minimal model under B — L, Z\ and Z% 



singlet considered in Ref. [8] is not needed for a scotogenic Majorana neutrino mass. 
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Figure 1: One-loop generation of Dirac neutrino mass in the minimal model. 



Whereas a scalar singlet was discussed as dark matter by itself many years ago p2J H7J HE] , 
our proposal may be considered a natural justification of its existence. 

Let the Yukawa interactions be given by fakVaN^rf and hkpNkvfjX ■ Without loss of 
generality, the A parameter of the trilinear Ax4>°ri term may always be chosen real, as well 
as the vacuum expectation value (<p°) = v. Let if = (tjr + ir]i)/y2, then there is a mixing 
between rjn and x°> but not between rjj and x°- Assuming in addition that rjj is a mass 
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eigenstate and denoting the mass eigenstates of the Vr) sector as £1,2 with mixing angle 
9, the one-loop Dirac neutrino mass matrix is then given by 

sin 9 cos 9 



{M v )„t 



16tt 2 v_ k 



9 2 2 2 

mi mi mi mi 

— ^ In == Ca s In 



m Ci - m k m N k m \ ~ m N k m N k 



(1) 



This is in complete analogy to that of the radiative Majorana seesaw [2], with suppression 
of the neutrino mass from the usual assumption of very large (now Dirac) as well as the 
loop factor. In addition, this diagram is only nonzero because of the soft breaking of Z^ . 
Thus, it is natural for the parameter A to be small. In the limit A = 0, the mixing angle 9 
in the above equation would be zero. 

We assume that there are three copies of (N, N c ) so that all three neutrinos obtain 
scotogenic Dirac masses. If there is only one copy, then two neutrinos will be massless, 
which is clearly unrealistic. If there are two copies, one will be massless, which is acceptable 
as far as present neutrino phenomenology is concerned. From Table 1, it can be easily 
confirmed that with three copies of z/ c , U(1)b-l will be anomaly-free. 

In this model, <3> is the SM Higgs doublet with the usual Higgs boson H as its only 
physical degree of freedom. It has the usual SM decay modes, except for corrections due to 
its interactions with rj and x°- F° r example, H may decay into C1C1 if kinematically allowed. 
If Ci is dark matter, this decay would then be invisible. It would affect the search for the 
SM Higgs boson, as studied already in Ref. [19]. Another possible effect is that the coupling 
of H to 7] + ri~ would change the one-loop decay of H to 77, thus affecting also the search for 
the SM Higgs boson via this channel. A third effect is the existence of the quartic XX^ ^ 
coupling, which may contribute significantly to the effective potential of H and modify its 
stability condition as a function of mass. It may also induce a one-loop contribution to 
the H 3 T term at finite temperature to cause a first-order phase transition needed for the 
electroweak baryogensis. 
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The couplings f a kL a Nj:ri contribute to radiative lepton flavor violating rare decays: 



m 



" 167T 

where 



with t = (m 2 Nk /rn^ + ). For t — > 0, t — > oo and i — > 1, the combination in the last parenthesis 
of Eq. ([3]) converges respectively to 1/6, l/(12t) and 1/24. For mjv fe 3> which is the 
seesaw limit, we find 

(E^f) 1/2 ~8xl0-(5<t±M) 1/4 GeV-. (4) 

We will consider first this scenario, so that the dark-matter candidate of our model is the 
lightest of the three exotic neutral scalars: £i,2 or rjj. 

The most general scalar potential consisting of $, r], and x is given by 

+ A 3 ($ t $)(^) + A 4 ($ t r7)(r ? t $) + ^A 5 ($^) 2 + H.c. 

+ \Kx^ + ^X 7 (&<5>) X 2 + -X 8 (rfv)x 2 + A X &V + H.c. (5) 

This potential preserves and breaks softly by the last term. The parameter A may 
be chosen real by a phase rotation of rj relative to $, but then A5 is in general complex. For 
simplicity, we choose it to be real so that r/j is a mass eigenstate and decouples from the 
(x°,Vr) sector. The resulting mass spectrum is given by 

m\ = 2\ lV \ (6) 
m l+ = vl + hv 2 , (7) 

™> 2 VI = I2 2 2 + (X 3 + X 4 -X 5 )v 2 , (8) 
(fij + Xrv 2 V2Av \ 



t In t t 2 - 5t 



2(t - l) 4 12(t — 1) £ 



(3) 



V2Av /4 + (A 3 + A 4 + X 5 )v 



This is very similar to previous studies such as Ref. [5] with an important conceptual differ- 
ence. Since the parameter A breaks z\ , it may be argued that it is small. This suppresses 
the radiative neutrino Dirac mass as well as the mixing between tjr and \. Hence the 
dark-matter candidate of this model can be dominantly a singlet and as a result, it can 
naturally evade the constraints from the electroweak interactions of a doublet. If there is 
no symmetry, the mixing between tjr and \ is then arbitrary, as in previous studies. 
Another difference is that r\i is not involved in the one-loop neutrino mass, contrary to the 
original Majorana case of Ref. [2]. The above possibility has also been discussed in [15]. In 
the following, we introduce a new possibility for dark matter candidate within the present 
scenario. 

Since m^ k are assumed to be very large in this scenario, the annihilation of the dark- 
matter scalars in this model do not proceed via their Yukawa interactions, but rather through 
their gauge or scalar interactions. Examples of the latter have been discussed extensively in 
the literature [IHl HZl UHl EOl EH [22j [23l El]. 

We consider next the lightest (call it N\) as dark matter. As shown previously [25J, 
this is subject to severe phenomenological constraints in the original model [2] of scotogenic 
Majorana neutrino mass. The reason is as follows. In order for N\N\ annihilation to account 
for the correct relic abundance, the rj masses cannot be too heavy and the Yukawa couplings 
f a k cannot be too small. However, these values are severely constrained by experimental 
upper limits on the \i — > rate, as already discussed. It is thus not a viable option, without 
some detailed fine tuning of parameters. To retain iVi as a natural dark-matter candidate, 
new interactions involving Ni need to be postulated, such as a singlet scalar [26J. In our 
present model, the h^N^i/jX Yukawa couplings are exactly what are required. They are not 
constrained by flavor- changing charged-lepton radiative decays, so they can be large enough 
for a realistic NiNi annihilation cross section to account for the relic abundance of dark 
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matter in the Universe today. In this scenario, the f ak Yukawa couplings as well as the A 
parameter are small and the mass of (j (which is mostly composed of x) is n °t much greater 
than m Nl . 

Combining (is, is c ) and (N, N c ) to form the four-component Dirac fermions v and N, their 
Yukawa interactions are given by 

£y = f ak Nk ( (u a rf - l a r] + ) + h kp N k CA? 1 ) vpX° + H.c. (10) 



where in this four- component notation, [(1 + 75)/2]i/ represents is c going backwards. For the 
dark-matter candidate N\, we assume h\p to be dominant, then 

aiNl + Nl _> Va + ^ = g _iL _ _ l_ < g (11) 

where to reach the last inequality we have used > m^. Similarly, 

\h* h* I 2 m 2 \h* h* I 2 

1 32 ^ W 1 +^ 2 1 ) 2 ^128^^ 1 ; 

Setting the sum of the two annihilation cross sections times the relative velocity equal to 
one picobarn, we find 

™*i < fc I^LM'j (1.4 TeV). (13) 

For \h\ a \ < 1, we then obtain < 4.2 TeV. With such light N\, the seesaw mechanism is 
not very effective. The smallness of the neutrino masses can be justified by the smallness of 
the trilinear A term which softly breaks Z% , and the smallness of the / Yukawa couplings. 
If the h couplings were not available, the cross section must have then come from the / 
couplings, which are restricted by /i — > ej, so the annihilation cross section would in general 
be too small for Ni to be a viable dark-matter candidate. If the B — L symmetry is gauged, 
there should be another annihilation mode N + N — > Z' — > is + v,l + l,q + q. This cross 
section is given by |27J 

4 2 

a = 7^ — 5 ^T? ( 14 ) 

ixv re i(Am z Ni - m%,y 

9 



The present lower bound on mz> from the Large Hadron Collider (LHC) [28] is estimated 



to be about 2 TeV. For g z > = ^J(5/8)g Y = 0.28 (i.e. the SO (10) limit), m z < = 2 TeV, 
and av re i = 1 pb, we find = 900 GeV. In this case, N\Ni production from Z' decay 
at the LHC is possible, as studied previously [29], except that N\ is now dark matter. It 
may however be inferred from the increase of the Z' invisible width on top of the expected 
Z' — > vv mode. As N\ is otherwise very difficult to produce, the existence of Z' seems to 
be the only realistic chance for it to be observed at the LHC, but still only indirectly. If r) + 
is light enough, it can be produced at the LHC. The subsequent decay of r] + into Ni and a 
charged lepton is a possible signature, as discussed in Ref. [30] . 

As for direct detection of dark matter in underground experiments, if B — L is not gauged, 
then N\ has no interaction with nuclei. If B — L is gauged, then the elastic scattering of N\ 
with nuclei may proceed through Z' exchange. The cross section per nucleon is given by [27] 

4m P 9Z> /, KN 

oo = —■ (15) 

7r m z , 

For gz> = 0.28 and mz> = 2 TeV, this implies <jq = 1.7 x 10~ 7 pb, which exceeds the 
XENON100 bound [3XJ of about 7 x 10~ 8 pb for m Nl = 900 GeV. This means that in this 
case, gz'/mz' should be reduced by a factor of 1.25 or more. 

This minimal model is also very suitable for the implementation of the non-Abelian 
discrete v4 4 symmetry [32J to the neutrino mass matrix [33J. In the charged- lepton sector, 
let (z/j, li) ~ 3 under A±, and either If ~ 1, 1', 1" as in Ref. [32] or l\ ~ 3 as in Ref. [31], then 
with $ ~ 3 or 3 + 1, and A4 breaking to the residual symmetry Z3, the charged- lepton mass 
matrix is diagonalized by the well-known unitary matrix 

/ll 1 \ 

1 u uj 2 , (16) 
\1 u 2 UJ J 

where u = exp(2ni/3). In the neutrino sector, let uf ~ 3, r) ~ 1, and x ~ 1 + 3, with the 
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soft scalar trilinear x$ V terms to break A4, the neutrino mass matrix becomes 

(a f e\ 

M u = fad. (17) 
ye d a J 

If e = / = 0, then neutrino mixing is tribimaximal, i.e. sin 2 9 12 = 1/3, sin 2 6*23 = 1/2, 
#13 = 0. This was known to be a good approximation of the measured neutrino mixing 
angles. However, two recent experiments have measured #13 to be definitely nonzero, i.e. 



sin 2 29 13 = 0.092 ± 0.016(stat) ± 0.005(syst) 

from the Daya Bay Collaboration [35], and 

sin 2 29 13 = 0.113 ± 0.013(stat) ± 0.019(syst) 

from the RENO Collaboration [36] . In that case, e and / should be nonzero. Let 

e-f s e+f 



(18) 



(19) 



(20) 



dy/2 ' dy/2 ' 

The parameters a, d, e, / are complex, and for small e, /, the eigenvalues of M. v are a + d, 
a, and a — d. We can always choose a to be real, the phase of d is then determined by the 
absolute values of the three masses [37]. For the small values of e and /, we find 



#13 = ^ 



v/3' 



tan 2 # 



12 



(1 - V2Re5) 2 + 2(ImS) 2 
2 l(l + Re5/y/2) 2 + (Im5) 2 /2_ ' 



(21) 



Thus, a nonzero #13 and a value of tan 2 #12 smaller than 0.5 can be obtained. More precisely, 

the neutrino mass matrix in the tribimaximal basis is now of the form 

/ mi me \ / a + d 5d \ 



^(1,2,3) 



m 6 m 2 m 5 
\ m 5 m 3 / 



5d a ed 
\ ed a — d J 



(22) 



If 5 = e = 0, the tribimaximal mixing is then recovered. This differs from the originally 
proposed deviation [33] for A±, which was updated recently [38], i.e. 

(mi mA I a + d- (6 + c)/2 zV3/2(c - b) \ 

m 2 m 5 = a + b + c \[2e . (23) 

\m 4 m 5 m 3 J \ iy/3/2(c — b) y2e a — d—(b + c)/2J 



W(l,2,3) 
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Given that m 4 = in Eq. ( I22p . we obtain the approximate relationship 

sin 2 20 2 3~l-8[ife([4 3 )] 2 . (24) 

Using the experimental bound sin 2 29 2 3 > 0.92, we find \Re(U e %)\ < 0.1. If we take the central 
value of \Ues\ to be 0.16 (corresponding to sin 2 26>i 3 = 0.1), we then obtain | tan5cp| > 1.3 
in this model. Details are given elsewhere [39J. 

Below we also mention briefly how a supersymmetric model of scotogenic neutrino mass 
may be constructed. Consider the superfield content listed in Table 2. There are two 
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Table 2: Assignments of the particles of the supersymmetric model under B — L, Z\ and 

zi B \ 



one-loop diagrams contributing to the Dirac neutrino mass as shown in Fig. 2. Note that 
supersymmetry is broken by the soft scalar trilinear X20?02 an d bilinear iViV c terms. There 




Figure 2: One-loop generation of Dirac neutrino mass in the supersymmetric case. 

are now many particles of odd Z\ as well as superpartners of odd R parity. There are thus 
at least two dark- matter candidates [ID]. Obviously the details of the dark sector are much 
more complicated. We will not study them further in this paper. 

In conclusion, we have studied a minimal model of radiative Dirac neutrino mass induced 
by dark matter. In order for the scotogenic Dirac neutrino mass to occur in one loop, we need 
to introduce a scalar singlet x° which mixes with the neutral component of a new electroweak 
scalar doublet (r] + ,r]°). It is thus a good theoretical justification for the existence of x°- in 
addition to the possibility of direct production at the LHC, the presence of rj + can modify 
the Higgs decay mode to 77. As shown in [JT], if the A3 coupling in Eq. (jSJ) is negative, it 
can lead to the enhancement of Br (if — > 77) in conformity of the recent observation at the 
LHC [12]. Moreover, the quartic coupling of x° with Higgs can stabilize its potential against 
radiative corrections. 

This minimal model also requires three heavy neutral Dirac fermions iVj. Depending on 
the mass spectrum, the dark matter might be either the lightest Dirac fermion N\ or one of 
the neutral scalars; i.e. the imaginary component rjj of rf or a linear combination of the real 
component tjr of rf and x°- in the latter case, depending on the mixing between t)r and x°, 
which should be small because of the soft breaking of , the annihilation rate due to the 
electroweak interactions can be made equal to about 1 pb which is a value dictated by the 
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dark matter abundance in the thermal dark matter scenario. 

If N\ is the dark-matter candidate, its annihilation can proceed via its Yukawa coupling 
with the right-handed neutrinos and x°- This is a possibility that does not exist within the 
scotogenic Majorana neutrino mass model because in that case the bounds from the // — > e^y 
constraints restrict the annihilation cross section of the Aq pair below the required value. At 
the LHC, Aq can then be produced via the decay of r] + and r}~ along with a charged lepton 

The B — L symmetry used to maintain the conservation of lepton number can be gauged. 
In that case, the present LHC lower bound on mz> is about 2 TeV. The interaction with the 
Z' boson provides another route for the annihilation of the Aq pair as well as a portal for 
the interaction with quarks and hence direct detection. The bound from the XENON100 
experiment already constrains the parameter space. 

This minimal model is also suitable for implementing an A4 symmetry in such a way 
that nonzero #13 and large Sep may be obtained. We have also briefly mentioned how a 
supersymmetric extension can be constructed. 

The work of E.M. is supported in part by the U. S. Department of Energy under Grant 
No. DE-AC02-06CH11357. The authors thank Galileo Galilei Institute for Theoretical 
Physics for its hospitality. YF acknowledges partial support from the European Union FP7 
ITN INVISIBLES (Marie Curie Actions, PITN- GA-2011- 289442). She is also grateful to 
ICTP for partial financial support and hospitality. 

References 

[1] T. Yanagida, in Proc. of the Workshop on Unified Theories and Baryon Number in 
the Universe (KEK, Tsukuba, Japan, 1979), edited by O. Sawada and A. Sugamoto, 



14 



p. 95; M. Gell-Mann, P. Ramond, and R. Slansky, in Supergravity, edited by P. van 
Nieuwenhuizen and D. Z. Freedman (North- Holland, Amsterdam, 1979), p. 315; R. N. 
Mohapatra and G. Senjanovic, Phys. Rev. Lett. 44, 912 (1980); P. Minkowski, Phys. 
Lett. 67B, 421 (1977). 

[2] E. Ma, Phys. Rev. D73, 077301 (2006). 

[3] C. Boehm, Y. Farzan, T. Hambye, S. Palomares-Ruiz and S. Pascoli, Phys. Rev. D 77, 
043516 (2008) |hep-ph/0612228] . 



[4] N. G. Deshpande and E. Ma, Phys. Rev. D18, 2574 (1978). 

[5] R. Barbieri, L. J. Hall, and V. S. Rychkov, Phys. Rev. D74, 015007 (2006). 

[6] See for example L. Lopez Honorez and C. E. Yaguna, JCAP 1101, 002 (2011) and 
references therein. 

[7] P.-H. Gu and U. Sarkar, Phys. Rev. D78, 073012 (2008). 

[8] Y. Farzan, Phys. Rev. D80, 073009 (2009). 

[9] E. Ma and D. Suematsu, Mod. Phys. Lett. A24, 583 (2009). 

[10] E. Ma, Phys. Rev. D80, 013013 (2009). 

[11] M. K. Parida, Phys. Lett. B704, 206 (2011). 

[12] D. Suematsu and T. Toma, Nucl. Phys. B847, 567 (2011). 

[13] S. Kanemura, T. Nabeshima, and H. Sugiyama, Phys. Rev. D85, 033004 (2012). 

[14] S. Kanemura, T. Nabeshima, and H. Sugiyama, Phys. Lett. B703, 66 (2011). 

[15] P.-H. Gu and U. Sarkar, Phys. Rev. D77, 105031 (2008). 

15 



[16 

[1?; 

[18 

[19 
[20 
[21 

[22 

[23 

[24 

[25 
[26 

[27; 

[28 



V. Silveira and A. Zee, Phys. Lett. 161B, 136 (1985). 
J. McDonald, Phys. Rev. D50, 3637 (1994). 

C. P. Burgess, M. Pospelov and T. ter Veldhuis, Nucl. Phys. B 619, 709 (2001) 
|hep-ph/0011335 . 



Q.-H. Cao, E. Ma, and G. Rajasekaran, Phys. Rev. D76, 095011 (2007). 



B. Patt and F. Wilczek, hep-ph/0605188 



V. Barger, P. Langacker, M. McCaskey, M. J. Ramsey-Musolf and G. Shaughnessy, 
Phys. Rev. D 77, 035005 (2008) [arXiv:0706.43TTl [hep-ph]]. 

S. Andreas, T. Hambye and M. H. G. Tytgat, JCAP 0810, 034 (2008) [arXiv: 0808 .02551 
[hep-ph]]. 

S. Andreas, C. Arina, T. Hambye, F. -S. Ling and M. H. G. Tytgat, Phys. Rev. D 82, 
043522 (2010) [ arXiv:1003.2595l [hep-ph]]. 

Y. Farzan, S. Pascoli and M. A. Schmidt, JHEP 1010, 111 (2010) [arXiv:1005.5"323l 
[hep-ph]]. 

J. Kubo, E. Ma, and D. Suematsu, Phys. Lett. B642, 18 (2006). 

K. S. Babu and E. Ma, Int. J. Mod. Phys. A23, 1813 (2008). 

S. Khalil, H.-S. Lee, and E. Ma, Phys. Rev. D81, 051702(R) (2010). 

G. Aad et at [ATLAS Collaboration], Phys. Rev. Lett. 107, 272002 (2011) 



[ar Xiv: 1108. 15821 [hep-ex]]. 
[29] K. Huitu, S. Khalil, H. Okada, and S. K. Rai, Phys. Rev. Lett. 101, 181802 (2008) 

16 



[30] Y. Farzan and M. Hashemi, JHEP 1011, 029 (2010) [arXiv:1009.0829l [hep-ph]]. 

[31] E. Aprile et al, Phys. Rev. Lett. 107, 131302 (2011). 

[32] E. Ma and G. Rajasekaran, Phys. Rev. D64, 113012 (2001). 

[33] E. Ma, Phys. Rev. D70, 031901 (R) (2004). 

[34] E. Ma, Mod. Phys. Lett. A21, 2931 (2006). 

[35] Daya Bay Collaboration: F. P. An et al, Phys. Rev. Lett. 108, 171803 (2012) 
jarXiv: 1203. 1669 [hep-ex]]. 

[36] J. K. Ahn et al. [RENO Collaboration], Phys. Rev. Lett. 108, 191802 (2012) 
jarXiv: 1204.0626 [hep-ex] ] . 

[37] E. Ma, Phys. Rev. D72, 073301 (2005). 

[38] E. Ma and D. Wegman, Phys. Rev. Lett. 107, 061803 (2011). 
[39] H. Ishimori and E. Ma. larXiv:1205.0075l [hep-ph]. 

[40] Q.-H. Cao, E. Ma, J. Wudka, and C.-P. Yuan. larXiv:0711.3881l [hep-ph]. 

[41] M. Carena, I. Low and C. E. M. Wagner. larXiv: 1206. 10821 [hep-ph], 

[42] CMS collaboration, in preparation; ATLAS collaboration, in preparation. 



17 



